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Abstract 


Coupling of sloshing dynamics within a partially filled rotating dewar of 
superfluid helium II with spacecraft dynamics are investigated in response to the 
realistic environmental disturbance forces and torques acting on the spacecraft 
during normal operation. This study investigate: (a) the rotating bubble of 
superfluid helium II reacting to combined environmental disturbances, including 
gravity gradient, aerodynamic, and magnetic forces and torques; (b) 
characteristics of slosh reaction forces and torques coupling with spacecraft 
dynamics; (c) the contribution of slosh dynamics to over-all spacecraft dynamics; 
and (d) activating of attitude and translation control system. The numerical 
computation of sloshing dynamics is based on the rotational frame, while the 
spacecraft dynamics is associated with non-rotational frame. Results show that 
the contributions of spacecraft dynamics are driven by the environmental 
disturbances coupling with slosh dynamics. Without considering the effects of 
environmental disturbances-driven slosh dynamics acting on spacecraft coupling 
with the spacecraft dynamics may lead to the wrong results for the development 
of spacecraft system guidance and attitude control techniques. 
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I. Introduction 


In order to carry out scientific experiments, some experimental spacecraft 
use cryogenic cooling for observation instrumentation and telescope, 
superconducting sensors for gyro read-out and maintain very low temperatures near 
absolute zero for mechanical stability. The approaches to both cooling and 
control involve the use of helium II. In this study, coupling of spacecraft 
dynamics with the realistic environmental disturbance forces and torques acting 
on the spacecraft during normal operation, as well as spacecraft control to 
maintain the spacecraft normal operation are investigated. The spacecraft adopts 
the cooling and boil-off from the liquid helium dewar as a cryogen and propellant 
to maintain the cooling of instruments, attitude and guidance controls and drag- 
free operation. Potential fluid management problems may arise due to asymmetric 
distribution of liquid helium and vapor or to perturbations in the liquid-vapor 
interface. Basic understanding in the coupling of slosh dynamics with six 
degrees of freedom in spacecraft dynamics plays a significant role in the 
development of spacecraft guidance and attitude control systems. 

Liquid helium II at a temperature of 1.8 K is used as the coolant. With 
its superfluid behavior, there is very small temperature gradient in the liquid 
helium. In the negligibly small temperature dependence of surface tension and 
also the small temperature gradient along the liquid-vapor interface which drive 
Marangoni convection, the equilibrium shape of the interface is governed by a 
balance of capillary, centrifugal, gravitational and dynamical forces. 
Determination of liquid-vapor interface profiles based on computational 
experiments can uncover details of the flow which can not be easily visualized 
or measured experimentally in a microgravity environment. 

The instability of the liquid-vapor interface can be induced by the 
presence of longitudinal and lateral accelerations and torques. Thus, slosh 
waves are excited, producing high and low frequency oscillations in the liquid 
helium. The sources of the environmental disturbance forces and torques for the 
operational spacecraft include aerodynamic force and torques, magnetic torques, 
and gravity gradient forces and torques 1 ' 3 . A recent study 4 suggests that the 
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high frequency accelerations may be unimportant in comparison with the residual 
motions caused by low frequency accelerations. 

The time -dependent behavior of partially-filled rotating fluids in reduced 
gravity environments was simulated by numerically solving the Navier-Stokes 
equations subject to the initial and boundary conditions 5 '*. At the interface 
between the liquid and the vapor fluids, both the kinematic surface boundary 
condition and the interface stress conditions for components tangential and 
normal to the interface were applied 30 . The initial conditions were adopted from 
the steady-state formulations developed by Hung et al 7 . Some of the steady-state 
formulations of interface shapes were compared with the available experiments 
carried out by Leslie 8 in a free-falling aircraft (KC-135). The experiments 
carried out by Mason et al 9 showed that the classical fluid mechanics theory is 
applicable to cryogenic helium in large containers with sufficiently large 
velocities 10 * 15 . 

At temperatures close to absolute zero, quantum effects begin to be of 
importance in the properties of fluids. At a temperature of 2.17 K, liquid 
helium has a A-point (a second-order phase transition); at temperatures below 
this point, liquid helium (helium II) has a number of remarkable properties, the 
most important of which is superfluidity. This is the property of being able to 
flow without viscosity in narrow capillaries or gaps. At temperatures other than 
zero, helium II behaves as if it were a mixture of two different liquids. One 
of these is a superfluid and moves with zero viscosity along a solid surface. 
The other is a normal viscous fluid. The two motions occur without any transfer 
of momentum from one to another for velocities below a critical velocity 10 15 . For 
the components of normal and superfluid velocities above a critical velocity, the 
two fluids are coupled through snarling in a complex tangle 10 ' 11 . 

The key parameter of critical velocity to distinguish one-fluid from two- 
fluid models is a function of fluid temperature and container size. In other 
words, in considering the dynamical behavior of helium II in a large rotating 
cylinder, a mixture of the superfluid and the normal fluid without separation of 
the two fluids for fluid velocities greater than the critical velocity is 
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accounted for in the model computation. Density concentration of superfluid is 
a function of temperature, which is also true for the surface tension and viscous 
coefficient for helium II 10 " 17 . In this study, the theory of viscous Newtonian 
fluids is employed with transport coefficients being a function of 
temperature 14,13 . 

In order to carry out the study of transient phenomena of coupling between 
slosh reaction torques and spacecraft dynamics, (a) Slosh dynamics based on fluid 
dynamics formulation, and (b) Spacecraft dynamics and control processes based on 
translational and rotational formulations, have been numerically solved 
simultaneously. In the meanwhile, it is true that spacecraft dynamics driven by 
slosh dynamics can cause the effect to lead spacecraft deviating from normal 
operation 3 . The results of present study disclose following behavior of 
spacecraft dynamics: (a) how differences in the dynamical characteristics of 
environmental disturbances, such as aerodynamic, magnetic and gravity gradient 
forces and torques contribute to the excitation of slosh waves in liquid-vapor 
systems, (b) how slosh dynamics provide time-dependent reaction force, torque and 
mass center fluctuations acting on the spacecraft normal operational systems; and 
(c) how slosh dynamics contribute to over-all time-dependent force, torque and 
mass center fluctuations acting on spacecraft during its normal operation under 
the actuation of attitude and translation control processes. 

II. Mon-Inertial Frame Mathematical Formulation of Sloshing Dynamics 

Experiment made by Andronikashvili 10 " 13 for rotating helium II shows that it 
is necessary to exceed a critical velocity for the interaction between the normal 
and superfluid components to establish entire bucket in rotation 10 ’ 13 . For the 
rotating dewar along outer diameter of 1.56 m and inner diameter of 0.276 m, the 
critical velocities to assume the interaction between the normal and super-fluid 
components are 6.4xl0' 7 and 3.6xl0*‘ m/s, respectively 10 " 13 . With rotating speed 
of 0.1 rpm, the linear velocities along the outer and inner walls of rotating 
dewar are 8.17xl0' 3 and 1.45xl0' 3 m/s, respectively, which are at least several 
hundred times on the order of magnitudes greater than that of the corresponding 
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critical velocities to assure the interaction between the normal and superfluid 
components of helium II 10 ' 13 . Based on this illustration, the problem under 
consideration have the special features to warrant an adoption of viscous 
Newtonian fluids formulation in this study. These special features are as 
follows: (a) Fluid velocities are at least several hundred times greater than 
that of the corresponding critical velocities; (b) These high fluid velocities 
can produce great enough vortex lines to snarl in a complex tangle to assure an 
interaction between the normal and superfluid components 10 ’ 17 ; (c) These vortex 
lines snarling with a complex tangle between the normal and superfluid components 
warrant the adoption of Newtonian fluid model 10 ’ 17 . 

Consider a closed circular dewar partially filled with helium II while the 
ullage is filled with helium vapor. The whole fluid system is spinning in the 
axial direction z of cylindrical coordinates (r, Q f z) , with corresponding 
velocity components (u, v, w) . The governing equations for non-inertial frame 
spacecraft bound coordinates spinning along its z-axis has been given in our 
recent studies 20-21 . In this formulation, all the environmental disturbances 
excluding gravity gradient forces are included in the g-jitter accelerations. 
In other words, dynamical forces, such as gravity gradient, g-jitter, and angular 
accelerations, and centrifugal, Coriolis, surface tension, viscous forces, etc., 
are given explicitly in the mathematical formulations 18 ’ 21 . In the computation of 
sloshing reaction forces, moments, viscous stress and angular momentum acting on 
the dewar wall of the spacecraft, one must consider those forces and moments in 
the inertial frame rather than in the non-inertial frame 20 ’ 24 . 

For the purpose of solving sloshing dynamic problems of liquid systems in 
orbital spacecraft under a microgravity environment, one must solve the governing 
equations 21 ' 22 accompanied by a set of initial and boundary conditions 23,24 . A 
detailed illustration of these initial and boundary conditions concerning the 
sloshing dynamics of fluid systems in microgravity was precisely given in Hung 
and Pan 20 ' 21 ' 23 . The computational algorithm applicable to cryogenic fluid 
management under microgravity is also given earlier 21 ’ 24 . Summarized computational 
algorithm are illustrated in Figure 1. In this study, in order to show a 
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realistic example, a dewar with an outer radius of 0.78 m, and an inner radius 
of 0.138 m f top and bottom radius of 1.10 m and a height of 1.62 m has been used 
in the numerical simulation. The dewar tank is 80% filled with liquid helium and 
the ullage is filled with helium vapor (total fluid mass is 287.6 kg). The 
temperature of liquid helium is 1.8 K. In this study the following data were 
used: liquid helium density = 145.7 kg/m 3 , helium vapor density = 1.47 kg/m 3 , 
fluid pressure = 1.6625 Pa, surface tension coefficient at the interface between 
liquid helium and helium vapor * 0.0353 N/m, liquid helium viscosity coefficient 
= 9.61x10’* m 2 /s; and contact angle = 0°. The initial profiles of the 
liquid-vapor interface for the rotating dewar are determined from computations 
based on algorithms developed for the steady state formulation of microgravity 
fluid management 7 ' 23 ’ 23 . 

A staggered grid for the velocity components is used in this computer 
program. MAC (marker-and-cell ) method 26 of studying fluid flows along a free 
surface is adopted. VOF (volume of fluid) method is used to track the dynamic 
behavior of liquid-vapor interface by solving finite difference equations 
numerically. Detailed computational algorithm are illustrated in Figure 1. The 
approximate flow velocity is calculated from the explicit approximation of 
momentum equations based on the results from the previous time step. Computation 
of pressure and velocity at the new time step are, thus, obtained from 
iteratively solving the pressure equation through conjugate residual technique 27 ’ 
30 . The configuration of liquid-vapor interface adjusted by the surface tension 
effect at the new time step are then finally obtained. The time step during this 
computation is automatically adjusted through the fulfillment of the stability 
criteria of computed grid size. Convergence criterion of the iteration of 
pressure equation is based on the computed velocity at each cell which satisfy 
continuity equation with the errors no more than 10* 3 of the velocity 
difference 31 . As for the volume conservation of liquid, a deviation of less than 
1 % error of volume is always guaranteed before a move to the next time step. 
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III. Characteristics of Environmental Disturbance Forces and Torques Acting on 

the Fluid Systems and Spacecraft 

To give an example of the GP-B Spacecraft, which is an Earth satellite 
orbiting at 650 km altitude directly over the poles, the orbit period, r 0 can be 
computed from following expression: 

j ? 3/2 

z o =2n —^Tn (*> ( 3 " 1 ) 

R e 9o 

where R, denotes the radius of Earth (= 6373 km); R c , the radius of the circular 
orbit (- R b + h = 7023 km); h, orbit altitude (= 650 km); and g Q , Earth gravity 
acceleration (= 9.81 M /s 2 ). For the case of GP-B, the orbit period r 0 = 97.6 min 
(5856 s), and orbit rate n = 2rr/r 0 = 1.07 x 10 -3 rad/s. 

(A) Gravity Gradient Forces and Torques 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 18 ‘ 21 

a ?g =-n 2 [3(r c d)r c -d] (m/s 2 ) (3-2) 

where A denotes gravity gradient acceleration vector; d, the vector (not an 
unit vector) from the fluid element to the spacecraft mass center; t c , an unit 
vector from the spacecraft mass center to the center of the Earth; and n, the 
orbit rate. 

To give an example of the GP-B Spacecraft, it is assumed that the gravity 
gradient exerted on the mass center of the spacecraft orbiting around the Earth 
on its specified orbit is zero. In other words, all the gravity acceleration 
exerted on the spacecraft is nothing but the gravity gradient acceleration which 
is defined in Equation (3-2). In this study, we are interested in investigating 
how gravity gradient acceleration affects the dynamical behavior of cryogenic 
fluid elements of helium and excite capillary effect governed liquid-vapor-solid 
interface disturbances. 

For the convenience of mathematical calculation, let us describe all the 
parameters involved in Equation (3-2) in terms of Cartesian coordinates. In 
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order to match with the computer simulation, mathematical derivation are 
considered in the first quadrant. Figure 2(A) illustrates the geometrical 
relationship of the parameters shown in Equation (3-2). 

For the case of GP-B, spacecraft is rolling about its boresight axis in a 
polar, 650 km altitude Earth orbit with the spacecraft boresight pointing at 
Rigel. The vehicle is rolling at a constant angular velocity about its boresight 
axis, while the boresight remains inertially fixed on Rigel. At time t = 0, the 
rolling axis of the spacecraft is aligned with the radial direction of the Earth 
center to the spacecraft mass center. Azimuth angle of Earth toward the location 
of the spacecraft mass center, i|r B , can be computed from the orbit period obtained 
from Eq.(3-1) under the normal operation of spacecraft 

2-Z-t (3-3) 

Z o 

where t is the time measured from the instant when the direction of spacecraft 
spinning axis is aligned with the radial direction of spacecraft mass center to 
the center of the Earth. Some modification is required if coupling of slosh and 
orbital dynamics in which the results of the spacecraft orbit deviating from 
normal operation are considered. 

Gravity gradient acceleration indicates that acceleration acting on any 
fluid mass inside the container increases two units of acceleration per unit of 
distance measured from the mass center of the container (point O g at Figure 2A) 
to the location of the fluid mass parallel along the radial axis from the mass 
center of the container to the center of the Earth (parallel to unit vector r c 
shown in Figure 2A) while the acceleration acting on the fluid mass decreases one 
unit of acceleration per unit of the shortest distance measured from the location 
of the fluid mass to the radial axis along the vector from the mass center of 
the container to the center of the Earth 1 *' 20 . 

On the other hand, gravity gradient torque m in the non-inertial frame 

g 

of spacecraft bound coordinates is given by 
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(3-4) 


M g = / ( a gg x d) dm = 3 n 2 [r c x ( Ir c ) ] ( N-m) 

where m (kg) denotes mass of the spacecraft dry mass, while J is the inertia 
tensor (kg.m 2 ) in the non-inert ial frame. For the case of GP-B with a mass of 
3000 kg, the dry mass inertia tensor in the body coordinate frame 32 is 

A [4777 0 O' 

J = 0 4777 0 (kg-m 2 ) (3-5) 

0 0 3 4 8 9 J 

Figure 3(A) shows the time variation of gravity gradient torque acting on the 
spacecraft systems for a full orbit period with component along (x,y,z) 
directions. It shows a harmonic oscillation in x- and y-directions in which 
peaks of the oscillations coincide with the quarter orbit periods of spacecraft. 
Gravity gradient torque acting on the spacecraft system along z-direction is 
relatively small and is on the 10‘ 2 order of magnitude smaller in comparison with 
that of the components in x- and y-directions. 

(B) Aerodynamic Forces and Torques 

Aerodynamic forces and torques are caused by molecular collisions of Earth 
atmospheric gases on exterior spacecraft surfaces, resulting in momentum and 
energy transfer. For the case of GP-B orbiting at 650 km above Earth in a 
circular, polar orbit with velocity vector magnitude v, the spacecraft encounters 
atmospheric gas with fluctuating density p. For a constant reference area A r and 
length L r , the aerodynamic forces and torques in the spacecraft bound frame are 
computed from 

F ai = jp v'A r C n (3-6) 

"ai = (3-7) 

where non-dimensional force C fi and torque coefficient vectors along direction 
i in the spacecraft bound frame represent the influence of spacecraft shape. 
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including aerodynamic shading effects , on the aerodynamic forces and torques, 
respectively. The force and torque coefficients are computed for all vehicle 
pitch and roll angles of the body frame. At 650 km altitude of free molecular 
flow region (mean free path - 2 X 10 5 m) , the momentum transfer is characterized 
entirely as interaction between incident gas molecules upon spacecraft surfaces 
which is a function of incident flow angle, molecular speed ratio, temperature, 
and surface roughness. The shading algorithm used in computing the aerodynamic 
force and torque coefficients sums elemental forces and torques only for those 
elements directly exposed to the incoming free molecular flow. Figures 2(C) and 
3(B) show time variation of aerodynamic forces and torques, respectively, along 
(x,y,z) directions with reference to Lockheed User's Handbook 32 extending in a 
full orbit period. It shows a series of harmonic oscillations with greater 
amplitudes of oscillations in terms of azimuth angle i|/ B in the ranges of -nl 4 < 
tJt, < rr/4, and 3rr/4 < t|/ B < 5rr/4 whereas smaller amplitudes of oscillations occur 
in the ranges of rr/4 < iJf B < 3rr/4, and 5tt/4 < 4f g < 7nr/4 for components along x- 
and y-directions . These results illustrate a fact that the spacecraft encounters 
a maximum drag when the incoming flow is along lateral or transverse direction 
of the spacecraft while the drag becomes minimum when the incoming flow is along 
the axial or longitudinal direction of the spacecraft for drag force and torque 
components along x-, and y-directions. As to the components of aerodynamic 
forces and torques along z-direction, it show a half cycle oscillation with 
respect to orbital period in which a negative value is resulted for the first 
half cycle (spacecraft translation along -z direction) whereas a positive value 
is shown for the second half cycle (spacecraft translation along +z direction) 
of the orbital period. 

(C) Magnetic Torques 

Magnetic torques on the spacecraft are due to interaction between the 
cryoperm shield and Earth magnetic field through the production of cryoperm 
magnetic dipole induced by the Earth magnetic field. Demagnetization 
coefficients characterize the vector relationship between the exterior magnetic 
field and the induced dipole for the shield. Although the cryoperm shield is 
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axially symmetric about the spacecraft boresight, unequal transverse axis 
demagnetization coefficients were used to construct the cryoperm shield dipole. 
An external magnetic torque M m (N*m) occurs when the GP-B cryoperm shield 
magnetic dipole vector D m (A*m 2 ) is in the presence of the Earth magnetic field 
B ( Tes la ) , namely , 

M m ^D m xB (N-m) (3-8) 

All vectors are in body frame coordinates* Figure 3(C) shows time variation of 
magnetic torques along (x,y,z) directions with referenced to Lockheed User's 
Handbook 32 extending in a full orbit period. It shows a series of harmonic 
oscillations while spacecraft completes a loop of orbit period along polar orbit 
across dipole geomagnetic fields of Earth along x- and y-directions . Magnetic 
torque acting on the spacecraft system along z-direction is relatively small and 
is on the order of 10" 2 order of magnitude smaller in comparison with that of the 
components in x- and y-directions. 

Characteristics of gravity gradient, aerodynamic and magnetic forces and 
torques acting on fluids and spacecraft systems are quite different. As to the 
dynamics of fluid systems, both gravity gradient and aerodynamic forces can 
activate fluid elements directly with the help of fluid dynamic equations, 
whereas the activation of gravity gradient, aerodynamic and magnetic torques on 
fluid elements have to be computed indirectly in fluid dynamic equations from 
angular velocities and angular accelerations which are coupled directly with 
spacecraft rotational equations. It is also interesting to note for differences 
between forces of direct contact, and forces of contact with distances (such as 
gravitational and magnetic forces) . In the present case of helium II fluids in 
rotating dewar, aerodynamic forces belong to forces of direct contact whereas 
gravity gradient forces are forces of contact with distances. There is no 
magnetic force acting on the helium II fluids directly because there is an 
interaction between cryoperm shield and geomagnetic fluid but not directly with 
helium II fluids. Comparison between the characteristics of gravity gradient and 
aerodynamic forces acting on the fluid masses show that gravity gradient forces 
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activate different magnitudes of forces on the fluid masses at different 
locations inside the dewar whereas aerodynamic forces drive same magnitude on 
fluid masses everywhere inside the dewar via forces of reaction through the 
container wall. 


IV Attitude and Translation Control Subsystem 
The following description of the GP-B Attitude and Translation Control 
Subsystem is a simplified version of LMSC/P08701, GP-B Attitude and Translation 
Control Subsystem Description 33 , provided by Lockheed in support of the GP-B 
design. 

(A) Attitude Control 

The attitude control system is based on a proven design used on the Hubble 
Space Telscope. The attitude control function generates the spacecraft torque 
command which is sent to the helium thruster control function. The nominal 
configuration uses a Proportional/Integral (PI) type control law with rate 
feedback and a roll rate filter which adds gain in the forward loop to attenuate 
disturbances at the roll frequency. Rate feedback is accomplished by using low- 
noise rate-integrating gyros and limiters are used to allow rate limited capture 
and maneuvers . 

(a) Inputs: 

co srr - Vehicle rate error (rad/s) 

0^ rr - Vehicle attitude error (rad) 

(b) Parameters: 

Control Gains 



Pitch/Yaw 

Roll 


Proportional Gain: 

Kp = 0.482 

Kp = 0.054 

(s' 3 ) 

Rate Gain: 

K r = 1.98 

K r « 7.26 

(s) 

Integral Gain: 

Ki = 0.25 

K t “ 0.0023 

(s' 1 ) 


Pitch/Yaw loop bandwidth » 0.15 cycles/s 
Roll loop bandwidth = 0.06 cycles/s 
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Constants are used to limit the value of some of the variables in the 
control processing logic. Limits of these constants are as follows: 


Pitch/Yaw 

(bw = 0.15 Hz) 

Roll 

( 

bw = 

0.06 Hz) 

Lml 

= 

1 

arcsec 

Liml 

= 

60 

arcsec 

Lm2 

_ 

0.01 

arcsec 

Lim2 

= 

0.05 

arcsec 

Lm3 

- 

1 

arcsec 

Lim3 


60 

arcsec 

Lm4 

- 

0.25 

arcsec 

Lim4 

= 

20 

arcsec 

Lm5 

= 

0.05 

arcsec 

Lim5 

= 

4 

arcsec 

Lm6 

3 X 

0.25 

arcsec 

Lim6 

= 

20 

arcsec 

Lm7 

3S 

0.25 

arcsec 

Lim7 

= 

20 

arcsec 


(c) Outputs 

Tcmd = Vehicle Torque Command 

(d) Processing: 

Calculations are applied to each axis (pitch, yaw, roll) every 0.1 s. The 
output is a vehicle torque command the designated axis updated at the 10 
Hz rate. In this calculation, x(n) refers to the value of x at the current 
control sample; and x(n-l) refers to the value of x at the previous control 
sample (delayed 0.1 s). Processing of calculations are as flows: 

1. Apply limit function to attitude error 

& mtt L - Attitude Error limited 
S«r_l = 

where Limit Function (Saturation) is: 

Liml(i) * i, if -Liml < i < Liml 
Liml(i) * Liml, if i > Liml; and 
Liml(i) = -Liml, if i < -Liml 

2. Compute integral path for attitude control 

6 Int (n) = Lim5[0 Int (n-l)] + K A 0 #rr x (n-l ) dt 

where 0 Int (n) - Integral path - attitude error for current time step; 
and dt * 0.1 s - control loop sample time 

3. Apply Roll Rate Filter to Attitude Control 
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(i) u r (n) - Roll filter input at the current time step 
u c (n) = 0 Int (n) + Lim4 [S. lr (n) ] + Lim6 [ K r &>. tr ( n ) ] 

(ii) x r (n) - Roll filter state at the current time step 
x r (n) = A r x t (n-1) + B t u t (n) 

(iii) y r (n) - Roll filter output at the current time step 
y r (n) = C r x r (n-l ) + D r u t (n) 

(iv) A r — Roll filter state transition matrix 



where filter coefficients: d x = -1.9999926, and d 2 = 0.9999937 

(v) B r - Roll filter input matrix 

fil 



(vi) C r - Roll filter output matrix 

C r~[^ n i~ n o d \) • ( n 2~ n o d 2)\ 

where filter coefficients: n 0 = 1.0031384, n 1 = -1.9999926, and 
n 2 * 0.9968553 

( vii ) Dr - Rollf ilter input / output coupling matrix 

M*o] 

4 . Compte Torque Command 

T^n) = K p ILim7[y r {n)\ 

where I - Vehicle Inertia about designated axis; and 

- Vehicle Torque Command for designated axis 
(B) Translational Dynamics 

The translational motion of proof mass with respect to the spacecraft (i.e. 
proof mass housing) is sensed by the capacitive pick-off sensor {drag-free 
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sensor) and used by the translational control function to control the spaceceaft 
translation. This is done to reduce inertial accelerations of the proof mass. The 
method demonstrates the method for determining the proof mass position error 
(i.e. - location of the proof mass housing with respect to the proof mass) . Steps 
of these methods are as folows: 

1. Sum all forces on the dry spacecraft (i.e. aerodynamics, magnetics. Helium 

s lo s h , • • . ) . 

2. Convert the total forces to inertial coordinates and solve translational 
equations of motion for the location of the dry spacecraft center of mass. 
The proof mass is assumed to be free from inertial forces and is located 
at the origin of the inertial coordinate frame. 

3. Compute the location of the proof mass housing relative to the proof mass 
by converting the center of mass to vehicle coordinates and adding the 
location of the proof mass housing relative to the spacecraft center of 
mass • 

^p«h_F» “ R<=»_p« 

where R^j* - location of the proof mass housing relative to the proof 
mass; ^ - location of the proof mass housing relative to the center of 
mass? and - location of the center of mass housing relative to the 

proof mass. 

(C) Trsanslational Control 

The translational control function generates the spacecraft force commands 
sent to the Helium thrusters. The normal configuration uses 
Proportional/Integral/Derivative (PID) type control law and a roll rate filter 
which adds gain in the forward loop to attenuate disturbances at the roll 
freqency. Positon feedback is accomplished by monitoring the proof mass positon 
with respect to t he spacecraft (proof mass housing) using a capacitive pick-off 
sensor. 

(a) Inputs: 

r vmh_pm ” location of the proof mass housing relative to the proof mass(m) 

(b) Parameters: 
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(i) Control Gains 

Each Axis (bw = 0.5 Hz) 

Proportional gain Kp = 2.44 
Rate gain k r - 1.2345 
Integral gain K L = 0.605 

(ii) Constants used to limit the value of some of the variables in the 
control processing logic. 

The limits are as follows: 

Each Axis (bw * 0.5 Hz) 

Liml = 5xl0“ 7 m 
Lim2 = 0.005 N 

(c) Outputs 

F^ - Vehicle Force Command 

(d) Processing: 

Calculations are applied to each axis every 0.1 s. The output is a vehicle 
force command for the designated axis updated at the 10Hz. rate. 

In this calculation: 

X(n) refers to the value of X at the current control sample; and 
X(n-l) refers to the value of X at the previous control sample (delayed 0.1 
s ) . 

Processing of calculations are as follows: 

1. Compute the proof mass error 

R mc - Proof mass position error 

R «rr = 

2. Compute integral path for translational control 

R Iat (n) = Liml [R Iot ( n-1 ) ] + K L R #xt (n-1 )dt 

where R Int (n) - Interal path - proof mass error for current time step; 
and dt = 0.1 sec - control loop sample time 

3. Apply Roll Rate Filter to Translational Control 

(i) u r (n) - Roll filter input at the current time step 
u r (n) = R Int ( n) + R.„(n) + K t [R. rr (n) - R. rr ( n-1 ) ] /dt 


15 



(ii) x r (n) - Roll filter state at the current time step 


x r (n) = A r x r (n-1) + B r u r (n) 

(iii) y r (n) - Roll filter output at the current time step 
y c (n) = C r x r (n-l ) + D r u r (n) 

(iv) A, - Roll filter state transiton matrix 


A r^ 


-<*1 - d 2 

1 0 


where filter coefficients: d x = -1.9999926, and d 2 = 0.9999937 

(v) B r - Roll filter input matrix 



(vi) C r - Roll filter output matrix 

C/\{n r n o d i)* ('* 2 -«o rf 2)] 

where filter coefficients: n„ = 1.0031384, n 2 = -1.9999926, and 
n 2 = 0.9968553 

(vii) Dr - Rollf ilter input / output coupling matrix 


M"o] 


4 . Compute force Command 

F^n) = Llm2[K p m-y r (n)] 

where m - Vehicle dry mass; and Fg^ — Vehicle Force Command for 
designated axis 


V Mathematical Formulation of Spacecraft Dynamics 
In spacecraft dynamics, a rigid body with six degrees of freedom, three 
being translational and three rotational, is considered. In this study, our 
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primary interest is to investigate the dynamical behavior of spacecraft driven 
by the coupling between each individual gravity gradient, aerodynamics, magnetic 
forces and torques, and the combined effects of these forces and torques coupling 
with spacecraft dynamics. 

The governing equations of three translational equations are given by 

— {mX,) = F Di + F li + F ci ( 5 - 1 ) 

where m, X, F D , F c and F L denote mass of spacecraft, inertial frame (non- 
rotational) coordinate, environmental forces such as gravity gradient and/or 
aerodynamic forces acting on the spacecraft, control force and slosh reaction 
forces (feedback from fluid system) acting on the spacecraft, respectively. 
Subscript i denotes component along direction i (=X,Y, or Z in non-rotational 
frame), while single dot and double dots at top (X r X) imply first and second 
order time derivatives of the parameter, respectively. Eulerian angles are 
defined to accommodate three rotational equations 34 . As usual, t|r, 6, and <p are 
defined as the heading, attitude and bank angles 34 . Three rotational equations 
in terms of Eulerian angles are given by 


J t ( 8+ijr 2 sin0cos0 ) +I a iircos0 (^-iJ/sinS) = M 0 

(5-2) 

D , * 

I a [ ( -0+{j/sin0+^&cos0 ) sin© -6 ( ^-ijrsin 0 ) cos© ] 
+J t (ilJcos 2 0 - 24 f 6 sin 0 cos 0 ) = M # 

(5-3) 

I a ( <p-tfsin9-i\i9cos9 ) = M ^ 

(5-4) 

Where 


M e = ( m lx +M dx +M cx) COS( P- (M L y +M Dy +M Cy ) sirup 

(5-5) 

= (M Lx +M Dx +M Cx )si.n<pcose+(M Ly +M Dy +M Cy )cos<pcose 
-(M Lz +M Dz +M Cz )sin8 

(5-6) 
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% = (* Lz + m Dz +m C z ) 


(5-7) 


where M Li , M oi and M ci are slosh reaction torque acting on the dewar and 
environmental disturbance torques (such as gravity gradient, aerodynamic and/or 
magnetic torques) and control torque acting on spacecraft, respectively, along 
direction i. M ti , M Di and are defined in non-inert ial (rotational) frame. 
With reference to Equation (3-5) for GP-B moments of inertia, I 1 = 1^, I 2 = I ^ 
and I 3 = I„. I* (= I 3 ) and I t (® I x = I a ) denote moment of inertial along axial 
and transverse directions, respectively. 

Activation of slosh dynamics in response to environmental disturbance 
forces and torques coupling with spacecraft dynamics certainly will induce 
angular velocity (co 1 ,c o 2 ,co 3 ) along (x,y,z) coordinates in rotational frame. Let 
( C0 1 / C0 2 f C0 3 ) = ( , & 2 r are tllQ t ^- me derivative of angular 

displacement along (x,y,z) coordinates, and qj denotes dewar rotating speed along 
z-axis. With the assumption of axially symmetric body about the z-axis, Eulerian 
equations shown in Equations (5-2) to (5-4) can be rewritten and expressed as 32 


= 7 - ("lx + M Dx + M cx ) 

0 2 +JC 2®l( d 3 + “) = -T^Ly + M Dy + M cy ) (5-8) 

1 2 

= y-(M Lz * m dz + M cz ) 

± 2 


in rotational frame, where K 3 = K 2 * (I l -I 3 )/I a , and K 3 =( Ij-Ii ) /I 3 . 

Both translational and rotational equations are initial-valued problems. 
Initial conditions for the translational equations are 


( X,Y,Z ) = (0,0,0) cm, and 
(X, Y t Z) = (0,0,0) cm/s at t = 0 


I 


(5-9) 


and that for the rotational equations are 
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(9,\\i f <p) = (0,0,0) rad, and 
(0,4r,0) = (0,0, 0.1) rpm at t = 0 


(5-10) 


i 


In this paper, rotational frame is adopted in the study of sloshing dynamics 
while non-rotational frame is chosen in the investigation of spacecraft dynamics. 
Conversion matrix A shall be used to make the conversion from the parameters in 
non-rotational frame (X,Y,Z) to that of the parameters in rotational frame 
(x,y,z) through the execution of following mathematical expression: 


X 

y 

z 



where 


A 


A A 
^xx ^XY 

Ay X Ayy 

A A 
^zx **ZY 



(5-11) 


(5-12) 


Ajm = -simlrcos# + cosijrsinSsin^, A** * sinilrsin^ + cosi(/sin0cos0 , A 1X = cos\J/cos0, 
* cosijfcos^ + sinijrs in8s in<p , A^ = -cosijjsin# + sim|/sin8cos<£ , A, y = sin\]icos8. 
Ax, = cos0sin<£, Ay, = cos8cos<£, and A,, = -sin8. 

As A is an orthogonal matrix, A“ x = A T , Equation (5-11) becomes 


X 

Y 

Z 



(5-13) 


Equations (5-11 to (5-15) shall be solved simtaneously to determine the 
values of the parameters of interest. Detailed algorithm of solving these 
equations coupling with slosh dynamics of fluid equation are illustrated in 
Figure 4 . 


VI Slosh Dynamics of Spacecraft 

Slosh dynamics driven by environmental forces and torques coupling with 
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spacecraft dynamics can be computed from Sections II, III, IV and V 
simultaneously. Figure 5(A) shows grid generation in which three-dimensional 
dewar, r-z plane at Q = 0°, and r-6 plane at z = L/2 are illustrated in Figures 
5(A-a), 5(A-b) and 5(A-c), respectively, where L is the height of dewar. Figure 
5(B) shows initial profiles of liquid-vapor interface with liquid fill level of 
80% in which three-dimensional, r-z plane at 9 = 0°, and r-Q plane at z = 108 cm 
are illustrated in Figures 5(B-a), 5(B-b) and 5{B-c), respectively. Here, 
initial interface profiles are computed under the conditions of g = 10“ 7 g 0 , co = 
0.1 rpm and \Jr B = 0°. Characteristics of slosh dynamics driven by each individual 
and combined environmental disturbance forces and torques coupling with 
spacecraft dynamics are investigated. 

Combined environmental disturbances, including gravity gradient, 
aerodynamic, and magnetic and control forces and torques, acting on fluid systems 
of dewar coupling with spacecraft dynamics are considered 35 ' 36 . Figure 6 shows 
the time sequence evolution of the three dimensional dynamical behavior of 
liquid-vapor interface oscillations driven by these combined environmental 
disturbances. It shows a series of asymmetric oscillations excited along liquid- 
vapor interface driven by these combine environmental disturbances. 

Figure 7 shows sloshing affected time evolution of fluid mass center 
fluctuations (in rotational frame) driven by combined environmental disturbances 
coupling with spacecraft dynamics. The values of fluid mass center fluctuations 
and their maximum absolute values are Max (Ax^, Ay^, Az^) = (1.47, 0.84, 5.6) 
mm, and Max (|x M |, (y^), |z M |) = (0.87, 0.44, 3.3) mm. These results imply Az^ 
> Ax^ > Ay^., and \z u .\ > | | > | | - 

Over-all slosh reaction forces and torques driven by combined environmental 
forces and torques have been investigated. Figure 8(A) shows the computed time 
variation of the fluctuations of slosh reaction forces (in rotational frame) 
exerted on the dewar in response to the combined environmental forces, and 
torques coupling with spacecraft dynamics. This figure shows the following 
results: (a) The values of slosh reaction force fluctuations and their maximum 
absolute values are Max (AF Lx , AF Ly , AF L2 ) * (3.6, 2.6, 5.7) 10'* N, and Max (|F Lx |, 


20 



|F ty |, |F t ,D = (2.0, 1.5, 2.9) 10- 4 N, respectively, which imply AF Ll > AF lx > AF ly , 
and | F »| > |F Lx | > | F ly | . ( b) Comparison of the characteristics of fluid mass 
center fluctuations and slosh reaction forces driven by combined environmental 
disturbances show that amplitude and fluctuations along the axial direction are 
greater than that along transverse directions. (c) Figures 9(A-a), 9(A-b), and 
9(A-c) show power spectral density with peak values of frequencies for 
autospectral analysis of slosh reaction forces exerted on the dewar in response 
to combined environmental disturbances along x, y and z directions, respectively. 
Major peaks shown are similar to the results concluded earlier. 

Figure 8(B) shows time variations of slosh reaction torques (in rotating 
frame) exerted on the dewar in response to the combined environmental forces and 
torques coupling with spacecraft dynamics. The values of slosh reaction torque 
fluctuations and their maximum absolute values are Max (AM Lx , AM Ly/ AM Lx ) = (1.6, 
1.7, 10" 3 ) 10' 4 Nm, and Max (|M Lx |, |M Ly |, |M Lx |) = (0.8, 0.7, 0.6xlCT 3 ) 10‘ 4 N m, 
respectively. This figure shows the following results: (a) AM^ - AM Ly » AM Lx , 
and | M lx | - |M ty | » | M Lx | . This means that the major slosh reaction torques 
driven by fluid systems exerted on the dewar in rotational frame are in the 
transverse rather than in the longitudinal directions, (b) Figures 9(B-a), 9(B- 
b) and 9(B-c) show power spectral density with peak values of frequencies for 
autospectral analysis of slosh reaction torques exerted on the dewar in response 
to the combined environmental forces and torques along x, y, and z directions, 
respectively. Major peaks shown are similar to the results concluded earlier. 

VII Control Force and Torque 

Figure 10(A) and 10(B) show the components of control forces and torques, 
respectively, along X, Y and Z direction based on the formulations of control 
equations illustrated in Section IV. 

Comparison between Figures 2(B) and 10(A) shows both characteristics of 
force variations and amplitude fluctuations for control forces and aerodynamics 
forces are basically similar except that there is out-of-phase for 180°. In other 
words, these two forces are equal in magnitude and opposite in direction. This 
also explains that major forces which drive spacecraft deviating from normal 
orbit are contributed by aerodynamic forces. In the meanwhile, fluid sloshing 
forces acting on the spacecraft also drive spacecraft deviating from normal 
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operation. Components of high frequency fluctuations imposed on the time- 
dependent variations of control forces are the results of product contributed by 
fluid sloshing forces acting on the spacecraft . 

Computed absolute values of maximum amplitudes for control forces are 
(|F Cx |, | f*cy j 9 | F c*|) = (2.6, 3.6, 2.5) 10‘ 3 N, while the saturation value (Lim2) 

of control subsystem documentation projects are 5xl0' 3 N. In other words, the 
amplitudes of control forces are not exceeding the saturation value during the 
processes of normal spacecraft operation based on this computation. 

Main environmental torques acting on the spacecraft are contributed by 
gravity-gradient, magnetic and aerodynamic torques. Based on the order of 
disturbances, these torques are on the order of 10“ 3 N *m. On the other hands. 
Figure 8(B) shows that the order of magnitude for slosh reaction torque is on the 
order of 10" 3 N m. This explains that the major deviation of spacecraft control 
shall be concentrated on the correction of combined effect resulting from 
spacecraft torque disturbances due to angular deviation from normal operation. 

VIII Resultant Forces and Resultant Torques 
Based on definition, total forces and torques acting in the spacecraft are 
= F l1 + F Di + F ci , and M i * M Li + M Dl + , respectively. Computed total forces 

and torques are illustrated in Figure 11(A) and 11(B). For the convenience of 
comparison, both F L and are transferred to rotational coordinates. 

Figure 11(A) shows that the order of magnitude for the total resultant 
forces are on the order of lO -4 N and are fluctuating from positive to negative 
values around the surrounding of 0 value when t>200 s. Computed results indicate 
that the absolute values of maximum amplitude of total resultant forces are 
(|F X |, | F y | , | F, | ) - (5.97, 4.23, 22.35)1CT 3 N. The range of maximum fluctuation 

of total resultant forces are (AF X , AF y , AF t ) * (9.15, 6.35, 4.30) 10' 3 N. The main 
cause of these maximum value fluctuations are induced by the fluid sloshing 
forces-related fluctuations at the beginning of computer simulation. In this 
computation, it is assumed that liquid helium is in static condition inside the 
rotating dewar in rotational frame at t = 0. Fluid is subjected to various 
environment forces which drive the fluid to non-equilibrium status in addition 
to the nature of incompressibility that make the initiation of large amplitude 
slosh oscillations which are fedback to spacecraft at t>0 . Fluid slosh damp out 
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gradually , and become the status of normal fluctuations after t>200 s. 

Similarly, characteristics of total resultant torques demonstrate the 
similar manner as that of total resultant forces illustrated earlier* The order 
of magnitude of total resultant torques are on the order of 10" 4 N *m. The 
absolute value of the maximum resultant torque are (|M X |, |My|, |M X |) 55 (2.62, 
1.86, 1.18)10~ 3 N m. The range of maximum fluctuations of total resultant torque 
are (AM X , AMy, AM X ) = (4.58, 3.45, 2.29)10" 3 N m. 


IX Attitude Errors 

Figures 12(A) and 12(B) show spacecraft angular rate errors and attitude 
errors, respectively. In the spacecraft normal operation, characteristics of 
angular velocity variation in rotating frame is (<o x , co y , co x ) = (0, 0, 0.1) rpm. 
To illustrate the errors of angular velocity. Figure 12(A) shows the value of 
spacecraft angular velocity errors in terms of Aw, = co x - co l0 , where co x denotes 
the actual value of the component of spacecraft angular velocity along z-axis 
while o> x0 shows angular velocity of spacecraft at normal operation (i.e., a) t0 = 
0 . 1 rpm ) . 

Absolute values of maximum angular velocity deviating from normal operation 
are (|(*) x |, | co y | , |Ag> x |) = (9.8, 22.1, 2.8)10"* rad/s. The range of maximum 
fluctuations of angular velocity deviating from normal operation are (Aco x , Aoo y , 
A(Aco x ) ) = ( 15.4, 38.4, 6.25) 10"* rad/s. In fact, these maximum values of 
fluctuations occur at the initial stage of the simulation. As indicated earlier, 
fluid inside the rotating dewar starts to receive environmental disturbances 
suddenly at t>0, this sudden induction of environmental disturbances introduce 
large amplitude fluctuations on the slosh dynamics and spacecraft system. Figure 
12(B) shows the value of spacecraft rotating angle errors in terms of A(p = <p - 
w l0 t, where co l0 denotes time variation of spacecraft normal operation with 
rotating angle <p , and t is time. In this expression, <p is actual value of angular 
displacement. During the spacecraft normal operation, = Q « 0. Figure 12(B) 
also illustrates the time dependent variations of parameters \\f , 6 and A<£ 
deviating from normal operation. The absolute values of maximum attitude errors 
are ( | \\f | , \Q \ , (A0 [ ) = (5.4, 114.8, 0.6)1 0 -7 rad, while the range of maximum 
fluctuations of attitude errors are (Aijj, A0, A(A<f>)) = (7.14, 30, 11.18)1 0 -7 rad. 

In this computation, sensor sensitivity to measure angular rate errors and 
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attitude errors are not included in our model. In fact, it is assumed that the 
corresponding control forces and torques will activate the spacecraft control 
systems based on control functions so long as there are angle errors and attitude 
errors no matter how small the amplitudes they are. This explains the magnitude 
of angle rate errors and attitude errors are so small. Similar results are also 
obtained for translational errors. The magnitudes of these errors will be 
amplified if sensor sensitivity and/or measurement errors are estimated and are 
taken into consideration. 


X Translational Errors 

Figures 13(A), 13(B) and 13(C) show translational acceleration, velocity 
and displacement errors, respectively, in non-rotating frame, with the activation 
of control system under normal operation. These figures show that the range of 
fluctuations in X and Y coordinates for translational errors are rather greater 
at the initial stage of simulation. As indicated earlier, spacecraft was driven 
by negative direction aerodynamic force suddenly at t £ 0 s, even though control 
force was acting on the system with delay of 0.1 s. This reaction at initial 
stage induces corresponding changes in velocity and in acceleration which induce 
helium II slosh dynamics and produce fluid reaction forces and torques fed back 
to the spacecraft system. These large amplitude slosh oscillations gradually damp 
out and the system becomes low amplitude oscillations after t > 200 s. Under the 
control system described, translational acceleration, velocity and displacement 
errors are relatively small and below the range of the effectiveness of control. 

Absolute values of maximum spacecraft translational accelerations deviating 
from normal operation are (|a x |, |a Y |, ]a g |) = ( 1.03, 2.03, 0*78)1 0 -4 cm/s 2 , while 
the range of maximum fluctuations of spacecraft acceleration deviating from 
normal operation are (Aa x , Aa Y , Aa x ) = (1.69, 3.70, 1.43) 10~ 4 cm/s 2 . Absolute 
values of maximum spacecraft translational velocities deviating from normal 
operation are <|V X |, |V T |, | V, j ) * (11.71, 2.84, 2.05) 10“ 3 cm/ s , while the range 
of maximum fluctuations of spacecraft velocities deviating from normal operation 
are (AV X , AV Y , AV t ) = (2.63, 5.61, 3.52)10" 3 cm/s. Absolute values of maximum 
spacecraft translational displacement deviating from normal operation are ( | X j , 

| Y | , | Z | ) = (1.03, 2.03, 0.78)1 0" 5 cm, while the range of maximum fluctuations 
of spacecraft displacement deviating from normal operation are (AX, AY, AZ ) - 
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(1.69, 3.70, 1.43) 10' 3 cm. 

As indicated earlier, computed results of translational acceleration, 
velocity and displacement will be amplified if sensor sensitivity and/or 
measurement errors are estimated and taken into consideration. 

XI Conclusions and Discussion 

It is shown in this simulation that the value of maximum variations of 
spacecraft attitude are on the order of 10" 6 rad while spacecraft translational 
errors are on the order of 10~ 3 cm, with the activation of control system under 
normal operation. The level of spacecraft accelerations are on the order of 10' 4 
cm/s 2 which is also equivalent to 10" 7 g 0 (g 0 = 9.81 m/s 2 ). In this simulation, 
with the neglection of initial stage and considering only for t>200 s, the level 
of spacecraft acceleration will be reduced to 10"* cm/s 2 which is also equivalent 
to 10“* g 0 . 

In this simulation, it is shown that control forces and torques required 
to carry over the necessary control processes are below the projected saturation 
values of the design requirements. 

The simulated values of spacecraft angular rate errors and spacecraft 
attitude errors will be amplified if sensor sensitivity and measurement errors 
are estimated and are included in the model computation. 

In this simulation, it is assumed that the fluid system was suddenly driven 
by environmental disturbances which cause relatively severe fluctuations of 
environmental forces and torques. This also implies that any large amplitude 
fluctuations can introduce large slosh dynamics which may cause a large 
disturbances which will be, deviating spacecraft from normal operation. 

In this simulation, it is shown that under the considered environmental 
disturbances, the spacecraft is under the safe side to be considered within the 
accuracy of spacecraft control system. 
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Figure Captions 

Computational algorithm for sloshing dynamics. 

(A) Geometry of gravity gradient acceleration, (B) Time evolution of 
aerodynamic forces. 

Time evolution of environmental disturbance torques acting on 
spacecraft during a full orbital period. (A) Gravity gradient 
torques, { B ) Aerodynamic torques, and (C) Magnetic torques. 
Computational algorithm for solving slosh dynamics of fluid equations 
in conjunction with translational and rotational equations of orbital 
dynamics and in coupling with translation and attitude control 
equations to carry out control processes. 

(A) Grid generation, (A-a) three-dimensional dewar, (A-b) r-z plane 
at 0=0°, (A-c) r-0 plane at z=L/2 <L=height of dewar); (B) Initial 

prof ile of liquid-vapor interface, (B-a) three-dimensional interface, 
(B-b) r-z plane at 0=0°, (B-c) r-0 plane at z=108 cm. 

Time sequence evolution of three-dimensional liquid-vapor interface 
driven by combined environmental disturnbanes including gravity 
gradient forces and torques, aerodynamic forces and torques, and 
magnetic torques, coupling with spacecraft dynamics. 

Time evolution of fluid mass center fluctuations along (x,y,z) 
directions (in rotational frame) acting on dewar driven by the 
combination of gravity gradient forces and torques, aerodynamic 
forces and torques, and magnetic torques. 

(A) Time evolution of slosh reaction forces (in rotational frame) 
along (x,y,z) directions driven by combined gravity gradient forces 
and torques, aerodynamic forces and torques, and magnetic torques. 

(B) Time evolution of slosh reaction torques (in rotational frame) 
acting on dewar driven by the similar combined environmental forces 
and torques. 

Power spectral density of slosh reaction forces and torques acting 
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including gravity gradient forces and torques, aerodynamic forces and 
torques, and magnetic torques. (A) Slosh reaction forces, and (B) 
Slosh reaction torques. 


Figure 10 Time evolution of control forces (in non-rotational frame) acting on 
spacecraft based on the translation and attitude control equations. 
(A) Components of control forces, and (B) Conponents of control 
torques . 

Figure 11 Time evolution of total resultant forces and torques acting on the 
spacecraft with the activation of control system under normal 
operation. (A) Total resultant forces and (B) Total resultant 
torques . 

Figure 12 Time evolution of spacecraft angular rate errors and attitude errors 
(in rotating frame) with the activation of control system under 
normal operation. (A) Spacecraft angular rate errors and (B) 
Spacecraft attitude errors. 

Figure 13 Time evolution of translational acceleration, velocity and 
displacement errors with the activation of control system under 
normal operation. (A) Translational acceleration errors, (B) 
Translational velocity errors and (C) Translational displacement 
errors . 
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Computational Algorithm for Sloshing Dynamics 



Fig. l 











(A) Geometry of Gravity Gradient Acceleration 
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Computational Algorithm for Solving Slosh Dynamics of Fluid Equations in 
Conjunction with Translational and Rotational Equations of Orbital Dynamics 
and in Coupling with Translation and Attitude Control Equations to Carry out 

Control Processes 



Fig. 4 




















(b) r - z plane (c) r - 0 plane 
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Time Sequence Evolution of Three-dimensional Liquid-vapor Interface 


(A) t = 0 s (B) t = 800 s 



(E) t = 3200 s 


(F) t = 4000 s 


Time Evolution of Fluid Mass Center Fluctuations 
(in rotational frame) 


(A) x Direction 



(B) y Direction 



(C) z Direction 



Fig. 7 
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Time Evolution of Slosh Reaction Forces (in rotational frame) Driven by 
Environmental Forces and Torques Coupling with Spacecraft Dynamics 

(A) Slosh Reaction Forces (B) Slosh Reaction Torques 
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Fig. 8 





Power Spectral Density of Slosh Reaction Forces and Torques Acting on the Dewar 
Driven by Environmental Forces and Torques Coupling with Spacecraft Dynamics 


(A) Slosh Reaction Forces 
(a) x Direction 



(B) Slosh Reaction Torques 
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Time Evolution of Control Forces and Torques Acting on Spacecraft Based on 
the Translation and Attitude Control Equations 

(A) Control Forces (B) Control Torques. 
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Time Evolution of Total Resultant Forces and Torques Acting on Spacecraft 
with Activation of Control System under Normal Operation 


(A) Total Resultant Forces (B) Total Resultant Torques. 
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Fig. 11 
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Time Evolution of Spacecraft Angular Rate Errors and Attitude Errors (in 
Rotating Frame) with Activation of Control System under Normal Operation 

(A) Spacecraft Angular Rate Errors (B) Spacecraft Attitude Errors 
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Fig. 12 














